Enhanced strange baryon production in Au+Au collisions compared to p+p at y/s^= 

200 GeV 
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We report on the observed differences in production rates of strange and multi-strange baryons in 
Au+Au collisions at ^/s'^= 200 GeV compared to p+p interactions at the same energy. The strange 
baryon yields in Au+Au collisions, when scaled down by the number of participating nucleons, are 
enhanced relative to those measured in p+p reactions. The enhancement observed increases with 

the strangeness content of the baryon, and increases for all strange baryons with collision ccntrality. 
The enhancement is qualitatively similar to that observed at lower collision energy ,/s^=17.3 
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GeV. The previous observations are for tlie bulk production, while at intermediate pr , 1 < Pt < 4 
GeV/c, the strange baryons even exceed binary scaling from p+p yields. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Nq,12.40.Ee 



One of the aims of studying relativistic heavy ion col- 
hsions is to observe how matter behaves at extremes of 
temperature and/or density. The energy densities in the 
medium produced by these collisions are far from that 
of ground state nuclear matter. Ultimately we hope to 
determine if they are sufficiently high to create a system 
where the degrees of freedoms are those of quarks and 
gluons, a state called the Quark-Gluon Plasma (QGP). 
By comparing the particles produced in A+A to those 
from p+p collisions, in which a QGP phase is not ex- 
pected, we can gain insight into the properties of the 
medium. 

Strange particles are of particular interest since the 
initial strangeness content of the colliding nuclei is very 
small and there is no net strangeness. This means that all 
strange hadrons must be formed in the matter produced. 
Originally, it was proposed that strangeness production 
would be increased due to the formation of a QGP com- 
pared to that from a hadron gas 1\. This enhancement 
is due to the high production rate of gg — > ss in a QGP, 
a process absent in the hadronic state. The subsequent 
hadronization of these (anti)strange quarks results in a 
significant increase in strange particle production, thus 
signaling a plasma was formed. 

The concept of enhanced strangeness production in the 
QGP can be recast in the language of statistical mechan- 
ics. A Grand Canonical Ensemble limit is likely only to 
be reached in the high multiplicity of heavy ion reac- 
tions. If this is the case, any measured enhancement is 
really a phase space suppression in p+p reactions that 
is removed in the heavy ion case. This lack of available 
phase space in small systems, such as those from p+p col- 
lisions, requires a Canonical ensemble to be used which 
results in a suppression of strangeness production when 
scaled to the appropriate volume 0, Hi- However, there 
is no a priori method for directly calculating this volume 
and thus the authors make the simplest hypothesis and 
assume that the volume is linearly proportional to the 
number of collision participants, (Npart)- The degree of 
suppression increases with the strange quark content of 
the particle. For sufficiently large volumes, the system 
is thermalized, the phase space suppression effects disap- 
pear, and the yields scale linearly with the volume, i.e. 
(Npart)- Initial measurements from the SPS suggested 
such a linear (A^nart) scaling 0|. However, it is not ob- 
served at RHIC or in the more recent SPS results Q. 

While the observables mentioned above are sensitive 
to the bulk of the produced particles with momenta be- 
low 2 GeV/c, further important information can be ex- 
tracted from intermediate and high pT particles. At 
RHIC, hadrons are suppressed at intermediate to high px 



when compared to binary-scaled p+p data at the same 
energy Q. This effect is attributed to the energy loss of 
partons as they traverse the hot and dense medium pro- 
duced 8, 9] . Measurements using identified particles help 
shed light on the details of the energy loss mechanism. 

In this paper we present further analysis of the high 
statistics measurements from p+p and Au-f Au collisions 
at ^s„„=200 GeV for strange and multi-strange baryon 
production at mid-rapidity as reported by the STAR col- 
laboration at RHIC [a, [ll I . Details of the STAR experi- 
ment are in _10| . Specific details of the trigger and detec- 
tors used to collect the data reported here can be found 
in [1, [llj and references therein. The Au-|-Au event 
sample consisted of 1.5x10^ central collision triggers and 
1.6x10® minimum bias triggers. The p+p results are 
from 6x10® minimum bias events. Particle identifica- 
tion is via the reconstruction of the charged daughter 
decay particles in the Time Projection Chamber. The 
decay channels used are A ^ p + tt~ , S~ ^ A -I- 7r~ 
p + TT^ + TT^ and ^ A + — > p + tt^ + plus 
the charge conjugates for the anti-particle decays. 

After cuts, to reduce random combinatorics, parent 
particles were selected if the calculated invariant mass 
fell within 3ct around the peak after background subtrac- 
tion. The data were corrected, as a function of pT, for 
efficiency and detector acceptance. Monte-Carlo studies 
showed that the corrections were constant as a function 
of rapidity over the measured regions. Further details 
of these reconstruction and correction techniques can be 
found in 0, [lH and references therein. Several contri- 
butions to the systematic uncertainty of particle yields 
were studied: detector simulation and efficiency calcula- 
tions, inhomogeneities of the detector responses, pile-up 
effects and the extrapolation of the data fits to zero pr- 
In p+p collisions an additional normalization error due to 
varying beam luminosity and trigger efficiencies of ~ 4% 
is included. The A yields were corrected for feed-down 
from multi-strange baryons using the measured spectra, 
the correction was of the order of 15%. 

For each species, i, the yield enhancement, E(i), above 
that expected from (Npart) scaling was calculated using: 



Yield^'^{i)(N^^ 



parti 



Yidd^«{z)(N^ait) 



(1) 



Fig. [T] shows E(i) as a function of (Npart), the inclu- 
sive proton data illustrate the effects for non-strange 
baryons l^]. Mid-rapidity hyperon yields measured as 
a function of centrality in Au-f-Au 5] and p+p [lH col- 
lisions were used. The number of participants, (Npart), 
and the number of binary collisions, (Ni,i„), were esti- 
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mated via a Monte-Carlo Glauber calculation [l^, 14 1. 
Since the p+p data were recorded with a trigger that was 
only sensitive to the non-singly difFractive (NSD) part of 
the total inelastic cross-section, all p+p yields have been 
corrected by o-^SD/'^inei' (=30/42), to obtain the total 
invariant cross-sections. 




10^ %a 



FIG. 1: (color online) Mid-rapidity E(i) as a function of 
(Npart) for A, A (|y| < 1.0), H", H+, Q~+ Ti'^ (\y\ < 0.75) 
and inclusive p (jyj < 0.5). Boxes at unity show statis- 
tical and systematical uncertainties combined in the p+p 
(p+Be) data. Error bars on the data points represent those 
from the heavy-ions. The solid markers are for Au-f Au at 
=200 GeV and the open symbols for Pb-I-Pb (|y| < 
0.5) at ^Sjyjy = 17.3 GeV [j. The arrows on the right axes 
mark the predictions from a GC formalism model when vary- 
ing rfrom 165 MeV {E(E-)^W.7, E(A)^2.6) to 170 MeV 
{E(S~)=7.5, E(k)=2.2). The red arrows indicate the predic- 
tions for H and the black arrows those for A , see text for 
details [li|. 



It can be seen that there is an enhancement in the 
yields over that expected from {Npart) scaling for all the 
particles presented. Since the proton yields are not cor- 
rected for feed-down, which is predominantly from the 
A and S, the p measurement is actually a sum of the 
primary protons and those from secondary decays. The 
integrated A -I- S'^ over inclusive p ratio varies from 30 to 
40 % for the p+p and Au-|-Au collisions respectively. If 
only primary protons were measured E (proton) would be 
closer to unity. A hierarchy in the scale of enhancements, 
which grows with increased strangeness of the baryon, is 
observed. This trend is predicted by Grand Canonical 
(GC) ensemble approaches, as is the fact that the E(i) 
values for each baryon/ anti-baryon pair are similar in 
shape 0- The difference in the scale of the enhance- 
ments for baryon and anti-baryon, especially at the SPS, 
is due to the existence of a non-zero net-baryon num- 
ber. However, the ratio of E(anti-baryon) to E(baryon) 
varies as a function of (Npart) at the SPS, possibly sig- 



nifying different production/annhilation mechanisms for 
(anti)particles at the SPS compared to at RHIC. For in- 
stance the net-A yields at the SPS, can be successfully 
described via multiple interactions of the projectile nu- 
clei [l5|. This effect is expected to be less significant 
at RHIC. It is also interesting to note that the mea- 
sured enhancements for the A, anti(2), and fl at RHIC 
are the same, within errors, as those calculated from the 
mid-rapidity SPS data (open symbols in Fig. [1]) despite 
an order of magnitude increase in the collision energies. 
Theoretical predictions using the GC ensemble approach 
predict a significant decrease in all the (anti)baryon en- 
hancements with collision energy |^]. A GC model, with 
a chemical freeze-out temperature of T—165 MeV and a 
baryon chemical potential, Hb =29 MeV calculates en- 
hancements of E('E.~ )=1Q. 7 and E(A)=2.6 for the most 
central Au-I- Au events at ys^=200 GeV [l3| • These en- 
hancement calculations cannot consistently describe the 
(anti)S and the (anti)A enhancements. However, the 
scales of the enhancements are very sensitive to the as- 
sumed freeze-out temperature and if T=170 MeV is used 
E{E.-)=7.5 and E(K)=2.2. 

While the measured enhancements are approximately 
constant for the inclusive protons they are clearly not for 
the A, 5, and Vl\ this is again counter to theoretical ex- 
pectations where the dependence of the strange baryon 
yields is expected to be linear with {Npart) for (Npart) ^ 
20. One explanation for this deviation from the theory is 
that the volume responsible for strangeness production is 
not linearly proportional to the geometrical overlap re- 
gion, as assumed in the model. A model that gives a rea- 
sonable description of the magnitudes and shapes of the 
enhancements with respect to centrality is described in 
[l3| . This model allows for an over-saturation of strange 
quarks, which varies with centrality, and thus does not 
invoke chemical equilibration. 

Fig. [1] is an average measurement of the difference in 
production between nucleus-nucleus and nucleon-nucleon 
collisions. Since the px distributions of the particles are 
approximately exponential these results are dominated 
by the physics occurring at < 2 GeV/c. Differences 
in the pr distributions for p+p and Au-|-Au data are 
studied by calculating the nuclear modification factor: 



Raa(pt, 



(FN'^^{i)/dpTdy 
TAAd^<J^^{i)/dpTdy' 



(2) 



Taa = (Nb^r^)/<ygel■ Fig- ^ shows Raa for A and the 
sum S+S for 0-5% Au+Au collisions along with those 
for inclusive p-l-p measurements [18. 19]. 

A striking feature of Fig. [2] is that both the central 
(top panel) and peripheral (bottom panel) Haa distribu- 
tions for the A and S~-|-S reach maxima that are much 
greater than unity, a value that would signify binary col- 
lision scaling. In fact the peripheral collision Haa distri- 
butions for the hyperons. Fig. ^p, are of approximately 
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FIG. 2: (Color Online) i?AA from a) 0-5% and b) 60-80% 
central Au+Au events for p-|-p [isl. [l9l|. A and H^+H^. Er- 
rors shown are statistical plus systematic added in quadra- 
ture. The band at unity shows the systematical uncertainty 
on (Nbin)- The dashed line below unity shows the expected 
value of Raa should the yields scale with (Npart) and the 
band around it shows the systematic uncertainty on (Npart)- 



the same magnitude as the central Kaa data, Fig. [2^, at 
intermediate to high px- These results are in contrast to 
the earlier reported suppression of high pT hyperons ob- 
served via Rcp [sl. [Til, fic 



in Fig. H 



Rcp{pT,i) 



2lj , these data are reproduced 



[(f iV^^"* {i)/dpTdy/ ( A^fc7„"* ) ] 
WNP-"P^ {i)/dpTdy/ {NlZ'"'')] ' 



(3) 



Non-strange hadrons reveal a similar suppression when 
using p+p or peripheral Au-|-Au colhsions as a reference. 
For pt > 1-5 GeV/c, unidentified charged hadrons show a 
suppression of the Au+Au spectra Comparing Haa, 
Fig. El to Rcp, Fig. [3l shows that Raa(A) Raa{^) 
^ Raa{p) but that Rcp (A) « Rcp(S) « Rcp(p), es- 
pecially at intermediate to high pT- This is possibly 
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FIG. 3: (Color online) Rcp Au-I-Au events for p-|-p 0-12%/60- 
80% [3, A+A and E--hH+0-5%/60-80% H. Also shown 
as the dashed curve are the resuhs for h++h" for 0-5%/60- 
80% (20l |. Errors shown are statistical plus systematic added 
in quadrature. The band at unity shows the systematical 
uncertainty on {Ntin}- The dashed line below unity shows 
the expected value of Rcp should the yields scale with (Npart) 
and the band around it shows the systematic uncertainty on 

(Npart) ■ 



due to phase space effects in the p+p data extending 
to this intermediate px regime. It is surprising that this 
decreased production in p+p events, while predicted in 
the soft physics/thermal production regime, i.e. pT < '2 
GeV/c, extends out to, and even dominates in, this in- 
termediate Pt region. Fig. ^ suggests that this effect 
is strong out to p-r 3 GeV/c. The shapes of the Rcp 
distributions at intermediate to high pT are generally in- 
terpreted as the result of parton energy loss in the hot 
dense matter and quark coalescence during hadroniza- 
tion. A comparison of Fig. EJi. and Fig. Eb shows that 
the turnover points occur at approximately the same pT- 
These data suggest that an enhancement of strangeness 
production has already set in in peripheral Au+Au col- 
lisions. This behavior is similar to that observed for the 
total yields in Fig. [U and quantitatively consistent with 
expectations from canonical suppression in p+p. Some 
portion of the Kaa peak may be explained via the Cronin 
effect, the observed increase in intermediate pr spectra 
in p-A collisions [23]. However, the Cronin enhancement 
stays constant, or possibly increases, as a function of cen- 
trality [2^, and this is not seen in our data. Effects 
due to radial flow in the Au+Au data are significant at 
RHIC energies, even for the multi-strange baryons 2J], 
but flow dominates only at low pr- The shapes of the 
Kaa distributions below 1 GeV/c are markedly differ- 
ent. The peripheral collision data indicate approximate 
binary scaling of the baryon yields while the most central 
data fall beneath binary scaling but significantly above 
that suggesting participant scaling. This again indicates 
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FIG. 4: (Color online) Comparison of Raa from data to HI- 
JING [27[ and EPOS ^ for 0-5% and 60-80% most central 
Au-I-Au collisions, for p-|-p , A and H~-|-H^. Errors shown are 
statistical plus systematic added in quadrature. The bands 
at unity shows the systematical uncertainty on {Nun}- The 
bands below unity on the left of the graphs are centered at the 
expected value of Raa should the yields scale with (Npart) 
and the widths of the bands indicate the systematic uncer- 
tainty on {Npart}- 



that there are different constraints on baryon production 
when going from p+p to peripheral to central Au-|-Au 
collisions. 

Comparisons to dynamical models can be used to 
understand in more detail how the close-to-equilibrium 
strangeness production can be achieved and whether the 
same mechanisms affect strange particle production at in- 
termediate Pt- In the HIJING model |25|| the yields and 
qualitative features of the strange baryon Raa measure- 
ments (solid curves in Fig. 3]) can only be obtained when 
baryon junctions and color strings are included 2^ 27 1. 



EPOS calculations [H,!!^ (dashed curves in Fig.g]) pro- 
duce similarly large differences in the hyperon Haa and 
Rcp [11] to those measured at RHIC and also give a 
qualitatively reasonable representation of the shape of 
the data. EPOS describes particle production via a par- 
ton model where Au+Au collisions are represented as 
many binary interactions. Each binary interaction is de- 
scribed by a longitudinal color field which is expressed 
as a relativistic string, or parton ladder. At a very early 
proper time, before hadronization, the collision region is 
split into two environments, the core, in which the den- 
sity of strings is high, and the corona, which surrounds 
the core and has a low string density. Production from 
the corona is due to collisions of nucleons at the periph- 
ery of the nuclei and modeled via string fragmentation. 
Corona production is thus similar to that from p+p colli- 
sions. Meanwhile particle production from the core is ap- 
proximated via a simple statistical hadronization process, 
similar to that described in [s^l , a collective flow profile is 
then imposed upon these particles. The relative weight of 
core to corona production varies with both centrality and 
particle species with the core dominating production in 
the most central events. Strange baryons are dominated 
by core production even in peripheral events. 

In summary, we observe enhanced strange baryon mid- 
rapidity production in Au-|-Au collisions, especially in 
the more central events, when compared to the (Npart) 
scaled p+p data from the same energy. The measured 
yields fail to scale with (Npart) as predicted if the GC 
regime is reached and the particle production volume 
scales with the geometrical overlap region. The magni- 
tudes of the suppressions are different to those predicted, 
but close to those measured at SPS energies. At interme- 
diate Pt the Raa values are higher than binary scaling 
of p+p data would predict. When attempting to under- 
stand the evolution of strange particle production from 
p+p to central Au-|-Au one must take into account both 
the effects due to a suppression of strangeness produc- 
tion in p+p and jet quenching plus quark recombination 
in Au-f Au collisions, with the former dominating at in- 
termediate PT- Since the measured Rcp values for all 
strange baryons equal those of the inclusive protons in the 
intermediate and high pT regions, and are significantly 
below binary scaling, the p+p-\ike suppression is already 
predominantly removed in peripheral Au+Au collisions. 
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